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The self-diffusion of two-dimensional small Ag islands (containing up to 10 atoms) on Ag(lll) 
surface has been studied using and self-learning kinetic Monte Carlo [J. Phys.: Condens. Matter 24, 
354004 (2012)] simulations. A variety of concerted, multi-atom and single-atom processes were au¬ 
tomatically revealed in these simulations. The size dependence of the diffusion coefficients, effective 
energy barriers as well as key diffusion processes responsible for island diffusion are reported. In 
addition, we have compared activation barriers for concerted diffusion processes with those obtained 
from Density Functional Theory (DFT) calculations. 

PACS numbers: 68.35.Fx, 68.43.Jk,81.15.Aa,68.37.-d 


I. INTRODUCTION 

Detailed understanding of diffusion processes is essen¬ 
tial for various surface related phenomena such as het¬ 
erogenous catalysis, epitaxial growth, surface reconstruc¬ 
tion, phase transition, segregation, sintering and other 
topics in surface science 1,2 . Therefore a great deal of ef¬ 
fort has been devoted towards the Self-diffusion of single 
atom and clusters on metal surfaces 3 both theoretically 
and experimentally 4-18 . There has also been both ex¬ 
perimental and theoretical efforts to find the activation 
barriers and prefactors for adatom diffusion process 19-31 
on various surfaces. Although various theoretical inves¬ 
tigations have been performed for Ag/Ag(lll) system 
using embedded atom method (EAM) 22,32 , corrected ef¬ 
fective medium theory (CEM) 33 , effective medium theory 
(EMT) 24, 33,34 and ab initio calculations 35,36 , but to the 
best of our knowledge there has been no systematic ex¬ 
perimental or theoretical effort to understand diffusion 
mechanisms and calculate their activation barriers for 
small Ag islands on Ag(lll) surface. Accordingly, in 
this article we focus on the identification and calculation 
of the diffusion processes and their activation barriers 
for self-diffusion of small Ag islands (1 — 10 atoms) on 
Ag(lll) surface as function of island size. 

Due to the presence of two three-fold adsorption sites 
(fee and hep) on fcc(lll) surface (Fig. 1(b)), adsorption 
energies for a single atom or cluster on this surface may 
differ for different systems leading to different diffusion 
pattern on this surface. It is therefore important to iden¬ 
tify these differences in adsorption energies and various 
activation energies for the diffusion processes in these sys¬ 
tems as accuratly as possible. Density functional theory 
(DFT) based calculations for the activation energies of 
diffusion processes, are most accurate to resolve such dif¬ 
ferences but are very time consuming. On the other hand 
classical molecular static calculations using interaction 
potential are very fast and easy to manage but their ac¬ 
curacy is determined by the interaction potential used. 
Therefore in this study we have done detailed DFT calcu¬ 
lations for the activation energies for some of the impor¬ 


tant diffusion processes and compare them with those ob¬ 
tained from interaction-potential based molecular static 
calculations to validate the accuracy of our approach. 

As potential energy surface of the atomically flat 
fcc(lll) surface is least corrugated, it results in low acti¬ 
vation barriers even for clusters to diffuse as a whole on 
this surface. Therefore studies of diffusion processes on 
this surface is a challenging problem for both experiments 
and simulations. For smaller islands (upto tetramer) all 
possible diffusion processes can be listed, but as island 
size gets larger it becomes more difficult to enumerate 
all possible diffusion processes a priori. On the other 
hand, as diffusion processes are rare events, molecular 
dynamics (MD) simulation cannot capture every micro¬ 
scopic processes possible. Therefore, to do a systematic 
study of small Ag island diffusion on Ag (111) surface 
we have used on-lattice self-learning kinetic Monte Carlo 
(SLKMC-II) method 37 , which automatically finds all the 
relevant atomistic processes and their activation barriers 
on-the-fly for system under study. 

This paper is organized as follows. In section II we 
briefly discuss SLKMC-II method and describe the way 
through which various diffusion processes are identified 
and their activation barriers are calculated. In section III 
we present details of single-atom, multi-atom and con¬ 
certed diffusion processes responsible for the diffusion of 
Ag islands of size 1 — 10 atoms on Ag(lll) surface. Sec¬ 
tion IV provide quantitative analysis of the results of the 
diffusion coefficients and effective energy barriers of the 
islands studied. Last section V is devoted to the discus¬ 
sion and conclusions. 


II. SIMULATION DETAILS 

To study self diffusion of Ag islands on fee Ag(lll) 
surface, we have carried out SLKMC simulations using 
an extended pattern recognition scheme 3 ' which includes 
both fee and hep sites in the identification of neighbor¬ 
hood around an atom. Further details on simulation cell 
size, calculation of energy barriers for various diffusion 
processes and other details can be found elsewhere 37 . 
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In order to check the accuracy of our approach for the 
calculation of energy barriers, we have verified the en¬ 
ergy barriers of some of the key concerted diffusion pro¬ 
cesses for various island sizes of Ag/Ag(lll) system using 
climbing image nudged elastic band (CI-NEB) method 
as implemented in the VASP code 38,39 employing the 
projector augmented wave (PAW) 39,40 and plane-wave 
basis set methods, setting the kinetic energy cutoff for 
plane-wave expansion to 500 eV and describing exchange- 
correlation interaction between electrons by the Perdew- 
Burke-Ernzerhof functional (PBE) 41 . For all cases, we 
have used slab size of 5x5x4 (for 10 atoms island we use 
a 7 x 7 x 4 slab) layers with a vaccume of 15 A. We re¬ 
lax all surface structures, using the conjugate-gradient 
algorithm 42 , until all force components acting on each 
atom is smaller than 0.01 eV/A. In all calculations, af¬ 
ter initial relaxation of the slab where all atoms were 
relaxed, we fix the bottom three layers. In calculating 
energy barriers and searching transition states, we first 
use the Nudged Elastic Band (NEB) method 43 for prelim¬ 
inarily determining a minimum energy path, then apply 
a Climbing Image Nudged Elastic Band (CI-NEB) 44 cal¬ 
culation. We find that using NEB followed by CI-NEB 
calculations is efficient in searching transition states than 
performing NEB alone. Here in, we use the term “(CI- 
)NEB” to refer a combination of NEB and CI-NEB cal¬ 
culations. In all (CI)-NEB calculations, we used 7 images 
in each case. 

In all our calculations we use the same pre-exponential 
factor of 10 12 s _1 for all diffusion processes for our KMC 
simulations, this has been proven to be a good assump¬ 
tion for the systems like the present one 45,46 . 


III. RESULTS 
A. General Details 

In order to uniquely identify whether an adatom is on 
an fee or on an hep site on the fcc(lll) surface requires 
knowledge of at least 2 layers, namely the top layer of 
the substrate and the layer below. Therefore, the same 
convention as in our previous study (Ref. 37) is followed. 
That is, the center of an upward-pointing triangle formed 
by the substrate atoms is an fee site, and the center of 
a downward-pointing triangle is the hep site . On the 
fcc(lll) surface, all atoms in an island can sit on fee 
sites or on hep sites or on a combination of both sites, 
called fee-, hep- and mixed-island, respectively. Further¬ 
more, activation barriers of single-atom processes depend 
on various factors including the type of step-edge, namely 
A- or B- type, (see Fig. 1(b) in Ref. 37) along which atom 
diffuses. Depending on the type of material either the 
fcc-island or the hep-island or mixed-island can be ener¬ 
getically favorable, but in case of Ag(lll) surface it is 
always the fcc-island with few exceptions as. 

A compact adatom island (called island from hereon) 
on an fcc(lll) surface can move in the three directions 



FIG. 1. (fee and hep sites on an fcc(lll) surface, with corre¬ 
sponding directions for concerted diffusion processes. 


(movement of a compact island along other three direc¬ 
tions going over the top sites usually results in processes 
with large energy barriers which are not selected for the 
range of temperature used in these simulations and hence 
we do not include them in this study to speed up simula¬ 
tions) shown in Fig.l. Note that the numbering scheme 
for the directions open to an atom on an fee site is in¬ 
verse to that for those open to an atom on an hep site 
(see Fig.l). We follow the enumeration convention for 
directions distinguished in Fig.l throughout the article 
in tabulating activation barriers for various processes. 

Our SLKMC simulations automatically discovered var¬ 
ious types of concerted, multi- and single-atom diffusion 
processes for the island sizes studied. 

Concerted diffusion processes involve all atoms moving 
together from all-fcc sites to all-hep sites or vice-versa. A 
concerted process can be either a concerted translation 
of the entire island in one of the three directions shown in 
Fig. 1 or concerted rotation of the entire island around 
the axis through the center of mass of the island per¬ 
pendicular to the substrate surface, either clockwise or 
anti-clockwise. Concerted rotational processes may (for 
example in the case of dimer) or may not (for example in 
the case of trimer) produce any displacement in the cen¬ 
ter of mass of the island whereas concerted translation 
processes always produces the maximum displacement 
(a/2^2) in the center of mass of the island. Note that 
in our simulations, concerted rotation processes were ob¬ 
served for dimer and trimer. On the other hand concerted 
translation processes were observed for compact islands 
of all sizes. For a particular island size (usually compact 
islands), activation barriers for the concerted processes 
along the three directions can be different (depending 
upon the type of step-edge along the diffusion direction) 
except for highly symmetric geometries like trimers and 
heptamers. Also activation barriers for concerted pro¬ 
cesses for a compact fcc-island along different directions 
are different than for those of the compact hep-island of 
the same size. We will discuss concerted diffusion pro¬ 
cesses for various island sizes in this section. 
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Various types of multi-atom diffusion processes like, 
shearing, reptation and rotation were found mainly for 
island sizes from 5-10 atoms. In the case of shearing pro¬ 
cess, part of the island (usually more than one atoms) 
move from fee to fee (or from hep to hep sites). On 
the other hand, reptation 4 ' is a two-step shearing pro¬ 
cess that move the cluster from all-fcc to all-hep sites 
or the reverse: first, part of the island moves from fee 
to hep sites (hep to fee) creating a dislocation; then the 
rest of the island moves from fee to hep (hep to fee). 
Hence at the intermediate stage, the island has mixed 
fee-hep occupancy. We have also observed multi-atom 
rotation processes, mostly involving trimers. Multi-atom 
processes like shearing is mainly found for compact or 
almost compact islands where as for mixed-islands, rep¬ 
etition and rotation were observed from island size 5 to 
10 atoms. Activation barriers for multi-atom (including 
shearing, reptation and rotation) processes depend on 
the type of the island (fee- or hep- or mixed-island) as 
well as type of step-edge along which diffusion process 
takes place. As island size increases, number of rnulti- 
as well as single-atom processes increases. Multi-atom 
processes are discussed in detail when we take up islands 
of size 8-10 where as single atom processes are discussed 
in subsection. Ill B 10. 


B. Diffusion Processes 

In this section key diffusion processes of the various 
general types (concerted, multi-atom and single-atom) 
are discussed in detail. In addition, activation barriers 
for some of the concerted diffusion processes are com¬ 
pared with those obtained from DFT calculations. Note 
that values on top of each pair of figures are the relative 
energies for those configurations. 


1. Dimer 

On fcc(lll) surface, a dimer can have three possi¬ 
ble arrangements; two dimers with both atoms on hep 
sites (called hh-dimer, Fig. 2(a)) or on fee sites (called 
ff-dimer, Fig. 2(b)) and a dimer with one atom on fee 
and the other on hep site (called fh-dimer, Fig. 3(a)). 
Note that both ff- and hh- dimers have equal lenght step 
edges of both A- and B-type. Our calculations show that 
ff-dimer is energetically more favorable (0.001 eV from 
EAM potential and 0.007 eV from DFT calculations) 
than hh-dimer whereas fh-dimer is less favorable (0.033 
eV) than ff-dimer. 

There are two different types of concerted diffusion 
processes, namely concerted translations (two symmet¬ 
ric translations along direction 2 and 3 and another one 
along direction 1) and concerted rotations (both clock¬ 
wise and anti-clockwise) for both hh- and ff-dimers as 
shown in Fig. 2 (a)-(d) and Fig. 2 (e) & (f), respectively. 
These processes transforms hh (ff) dimer into ff (hh) 






FIG. 2. Different configurations for dimer along with activa¬ 
tion barriers for concerted processes (values in brackets are 
from DFT calculations), (a & b) All atoms on hep sites (hh- 
dimer) and all atoms on fee sites (ff-dimer) showing concerted 
translation processes (c & d) Same as (a & b) but showing 
second type of concerted process, (e & f) Concerted rotation. 


dimer. In the case of hh-dimer, activation barrier for con¬ 
certed translation (Fig. 2(a)) along direction 2 (Fig. 1(b)) 
is 0.115 eV where as along direction 1 (Fig. 1(b)) it is 
0.097 eV (Fig. 2 (c)) while for concerted rotation around 
the top site ( Fig. 2 (e)) it is 0.055 eV. For the case 
of ff-dimer, activation barrier for concerted translation 
(Fig. 2(b)) along direction 2 (Fig. 1(b)) is 0.116 eV where 
as along direction 1 (Fig. 1(b)) it is 0.098 eV (Fig. 2 (c)) 
while for concerted rotation (along direction 2 - anti¬ 
clockwise) around the top site ( Fig. 2 (e)) it is 0.056 
eV. Note that, due to symmetry, activation barrier along 
direction 3 for concerted translation and rotation (clock¬ 
wise and anti-clockwise) is the same, for both hh- and 
ff-dimer. 

Our database also contains single atom processes which 
converts ff-dimer to fh-dimer (mixed-dimer) with energy 
barrier of 0.055 eV (Fig. 3(b)) with reverse barrier of 
0.021 eV (Fig. 3(a))to convert back to ff-dimer. Simi¬ 
larly an energy barrier of 0.055 eV converts an hh-dimer 
(Fig. 3(d)) to fh-dimer with a reverse barrier of 0.022 
eV (Fig. 3(c))to go back to hh-dimer. Our database also 
contains single atom long jump processes (in which an 
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TABLE I. Activation barriers (eV) of concerted dimer trans¬ 
lations processes. Values in brackets are from DFT. 


Directions 

fee 

hep 

1 

0.098 

0.097 

2 

0.116(0.127) 

0.115(0.120) 

3 

0.116(0.127) 

0.115(0.120) 





I 



FIG. 4. Two compact trimers obtained from fee dimer, (a) 
F3H, all atoms on fee sites centered around an hep site (b) 
F3T, all atoms on fee sites centered around a top site. 


FIG. 3. Different configurations for dimer along with activa¬ 
tion barriers for single atom processes, (a & b) Single atom 
process to transform mixed dimer (a) to ff-dimer and back, 
(c & d). Single atom process to transform mixed dimer (c) to 
hh-dimer and back, (e & f) Long jump single atom processes 
including detachment. 

atom of dimer moves from fcc(hcp) to next fcc(hcp) site) 
as well as detachment processes as shown in Fig. 3(e & 
f). We found that both ff- and hh-dimers diffuse via con¬ 
certed (both rotaion and translation) as well as single 
atom processes whereas fh-dimer diffuses via single atom 
processes. 


2. Trimer 

Two different types of compact triangular trimers can 
be obtained from each ff- and hh-dimers of Fig. 2(b & 
a). Two types of compact triangular fee timers can be 
obtained from an FF dimer by attaching a third atom 
to the sites mentioned by black open circles in Fig. 4(b): 
one centered around hep site (called F3H trimer with all 


step edges of B-type as shown in Fig. 4(a)) and the other 
centered around top site (called F3T trimer with all step 
edges of A-type as shown in Fig. 4(c)). Similarly two 
types of compact triangular hep trimers can be obtained 
from HH dimer: one centered around fee site (called H3F 
trimer with all step edges of B-type) and the other cen¬ 
tered around top site (called H3T trimer with all step 
edges of A-type) as shown in Fig. 5(c & b). 

In the case of compact F3T and H3T trimers two types 
of concerted processes were found, non-diffusive symmet¬ 
ric rotations (both clockwise and anti-clockwise as shown 
in Fig. 5 (e & f)) and diffusive translations (symmetric 
along all three directions) as shown in Fig. 5 (b & d) 
along direction 2 for an F3T trimer. Activation bar¬ 
rier for the concerted rotation process for F3T trimer 
is 0.149 eV (Fig. 5(e)) while for H3T trimer it is 0.147 
eV (Fig. 5(f)). Concerted rotation process transforms an 
H3T timer into an F3T trimer and vice versa. In case of 
a concerted translation process, the activation barrier is 
0.180 eV and 0.172 eV for F3T and H3T trimer respec¬ 
tively (Fig. 5(d & b)). Concerted translation processes 
transform an F3T(H3T) timer into an H3F(F3H) timer 
as shown in Fig. 5 (a - d). 

Different than the case with F3T and H3T Trimers, 
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FIG. 5. Arrangement of atoms in a compact trimer along 
with concerted processes and their activation barriers. Values 
in brackets are obtained from DFT based (CI)-NEB calcula¬ 
tions. (a) F3H, all atoms on fee sites centered around an hep 
site (b) H3T, all atoms on hep sites centered around a top 
site (c) H3F, all atoms on hep sites centered around an fee 
site (d) F3T, all atoms on fee sites centered around a top site 
(e) & (f) concerted rotational process of F3T and H3T and 
their activation barriers 


F3H and H3F have only concerted translation processes 
in all three directions with the same energy barriers 
(due to symmetry) of 0.179 eV and 0.171 eV respec¬ 
tively (Fig. 5(a & c)). These concerted diffusion pro¬ 
cesses transform F3H to a H3T timer and H3F to a F3T 
trimer. Energy barriers obtained from DFT calculations 
for the concerted translations for the case of F3H and 
H3T are 0.196 eV and 0.166 eV respectively which are in 
close agreement with the values obtained from our static 
calculations. Our results for the activation energies of 
various concerted processes agrees very well with those 
obtained from Effective Medium Theory (EMT) poten¬ 
tial as reported in 48 . 

Our simulations also finds single atom shape changing 
processes for trimers. For compact trimers (Fig. 6 (a-d)), 
an atom can can move in 4 different directions as shown 
in Fig. 6. Directions 1 & 4 correspond to edge diffusion 
process which opens up the trimer (into an non-compact 
trimer) while direction 2 & 3 correspond to detachment 
process. We note that these processes move atoms from 
fee (hep) to nearest fee (hep) site. Activation barriers 
in these 4 directions for different types of trimers are 
shown in Fig. 6. Other types of single atom processes 
like corner rounding to convert non-compact trimers to 
linear trimers (and back) are also automatically found 
in our simulations. We note that as the island size gets 
bigger, types of these single atom processes becomes large 
and we categorized them in single atom processes section 
for cluster sizes lager than 4 atoms. It can be seen from 
Fig. 6 that single atom processes for trimer have high 



FIG. 6. Single atom processes possible for a trimer and their 
activation barriers based on type of trimer along different di¬ 
rections. 



FIG. 7. Possible shapes of compact tetramers obtained from 
a compact F3H trimer. (b) Compact F3H tetramer - dotted 
circles show three possible positions marked as a, b and c for 
the fourth atom - to get three compact shaped fcc-tetramers. 
(a & c) Parallelogram shaped compact fcc-tetramers. (d) Di¬ 
amond shaped compact fcc-tetramer. 

activation barriers and were rarely observed during KMC 
simulations. 


3. Tetramer 

Adding another atom to three different available po¬ 
sitions to the compact F3H trimer shown in Fig. 7 (b) 
results in the formation of three compact fcc-tetramers 
with a long (along the line joining the farthest atoms) 
and a short diagonal (perpendicular to the long diagonal 
) as shown in Fig. 7 (a, c & d). Similar tetramers are ob¬ 
tained from F3T trimer. Note that both F3H and F3T 
trimers results in the same three compact fcc-tetramers ( 
Fig. 7 (a, c & d)) which are symmetric to eachother. Sim¬ 
ilarly both H3F and H3T trimers results in three compact 
hep-tetramers with geometry similar to those of compact 












6 





FIG. 8. Various diffusion processes for compact 
tetramers.Values in brackets are obtained from CI-(NEB) 
based DFT calculations, (a & b) Concerted diffusion along 
direction 2 for a compact fcc-tetramer (hcp-tetramer). (c 
& d) Multi-atom dimer shearing process for a compact fcc- 
teramer (hcp-tetramer). 


TABLE II. Activation barriers (eV) of concerted processes for 
compact fee- and hep-tetramers as shown in Fig. 8(a & b) 


Directions 

fee 

hep 

1 

0.193 

0.191 

2 

0.193 

0.191 

3 

0.245 

0.244 


fcc-tetramers shown in Fig. 7 (a, c & d). It should be 
noted that all compact fee- and hep-tetramers have 2 A- 
and B-steps of same length. 

Our calculations show that fcc-tetramer is 0.002 eV 
more stable that hcp-tetramer. Three types of concerted 
processes are possible for a compact tetramer; one along 
each direction as shown in (Fig. 1(b)). Two out of the 
three concerted processes (along direction 2 & 3) for 
each fee- and hcp-tetramer are symmetric processes. Our 
SLKMC simulations also found single atom as well as 
multi-atom processes for tetramers. 

Fig. 8(a & b) shows concerted diffusion process, which 
converts an fcc-tetramer to an hcp-tetramer, along di¬ 
rection 2 (barrier along direction 3 is also same due 
to symmetry) for an fcc-tetramer with energy barrier 
of 0.193 eV with a reverse energy barrier (for an hcp- 
tetramer to covert to fcc-tetramer) of 0.191 eV whereas 
energy barrier along direction 1 is 0.245 eV with reverse 
barrier (for hcp-tetramer) of 0.243 eV (see Table. II). 
Note that concerted processes in symmetric directions 
are not shown. We also found multi-atom process as 
shown in Fig. 8(c & d) in which two atoms (along one 
of the A-steps in the tetramer) move together in the 
same direction like a shearing mechanism. This diffu¬ 
sion process has an activation barrier of 0.285 eV and 
0.276 eV for the compact fee- and hep-tetramers respec¬ 



FIG. 9. Possible shapes of compact fcc-pentamers obtained 
from a compact fcc-tetramer (a) Compact-fcc tetramer - sites 
with dotted circles marked as b and c for the fifth atom - 
to get compact fcc-pentamers: (b) with three B-steps (one 
B-step edge of length 2 and two B-step edges each of length 
1) and one A-step (of length 1). (c) three A-step edges (one 
A-step edge of length 2 and two A-step edges each of length 
1) and one B-step edge (of length 1). 


tively. This multi-atom process converts parallelogram 
shaped fcc-tetramer (Fig. 8(c)) to diamond shaped fcc- 
tetramer (Fig. 8(d)) and vice versa. Anothe multi-atom 
process involving two atoms (along the other A-step of 
the tetramer in Fig. 8(c)), called dimer shearing, converts 
it into parallelogram shaped fcc-tetramer of Fig. 8(a). 
Similar multi-atom processes are also present for compact 
hcp-tetramer. These dimer shearing processes are also 
present in island sizes greater than tetramer. Note that 
these multi-atom processe are crucial for random motion 
of tetramer island, their absence results in anomolous dif¬ 
fusion of tetramer. Multi-atom shearing processes along 
B-steps of the fee- and hep-tetramers have very high en¬ 
ergy barries and hence we do not report them here. 


4. Pentamer 

Two compact fcc-pentamers can be obtained by at¬ 
taching an atom to the compact fcc-tetramer of Fig. 8 (a) 
as shown in Fig. 9(a). Although the geometry of com¬ 
pact fcc-pentamers thus obtained looks similar, attach¬ 
ing an atom to an A-type step edge of an fcc-tetramer 
results in an fcc-pentamer with three A-type steps (one 
A-type step of length 2 and two A-type steps each of 
length 1) and one B-type step (of length 1) as shown 
in Fig. 9(c) (called A-type fcc-pentamer from hereon) 
whereas attaching an atom to a B-type step edge of an 
fcc-tetramer results in an fcc-pentamer with three B-type 
steps (one B-type step of length 2 and two B-type steps 
each of length 1) and one A-type step (of length 1) as 
shown in Fig. 9(b)) (called B-type fcc-pentamer from 
hereon). Similarly two compact hep-pentamers called A- 
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type hcp-pentamer (with three A-type steps and one 13- 
type step) and B-type hcp-pentamer(with three B-types 
and one A-type step) can be obtained by attaching an 
atom either on A-step or B-step of the hcp-tetramer (see 
Fig. 10(b & d)). Compact A-type fcc-pentamer as shown 
in Fig. 10(c)) is 9 meV energetically more favorable than 
B type hcp-pentamer (Fig. 10(d))) where as B-type fcc- 
pentamer (see Fig. 10(a))) is 5 meV less favorable and 
A-type hcp-pentamer (shown in Fig. 10(b))). 

In our SLKMC simulation we found that compact A- 
and B-type fee- and hep-pentamers diffuses mostly via 
concerted diffusion processes which displace the island 
as a whole from fcc-to-fcc (hep-to-hep) and vice-versa. 
Fig. 10 shows an example of concerted diffusion process 
and corresponding activation barriers for parallelogram 
shaped compact A- and B-type hep- (Fig. 10(b & d)) 
and fcc-pentamers (Fig. 10(c & a)) along direction 3. Ta¬ 
ble III shows activation barriers for concerted processes 
for a compact parallelogram shaped A-type and B-type 
fee- and hep-pentamers (see Fig. 10 (a)&(b)) in all 3 
directions. As can be seen from Fig. 10 (a & b) that 
fcc-to-hcp concerted diffusion process for parallelogram 
shaped B-type fcc-pentamer have lower activation barri¬ 
ers than hep-to-fee diffusion process for compact paral¬ 
lelogram shaped A-type hcp-pentamer as compact par¬ 
allelogram shaped A-type hcp-pentamer is 5 mEV more 
stable than compact parallelogram shaped B-type fcc- 
pentamer. Similarly compact parallelogram shaped A- 
type fcc-pentamer is 8 meV more stable than compact 
parallelogram shaped B-type hep-pentmaer as shown in 
Fig. 10 (c & d). Also note that concerted processes in 
other two directions (not shown in Fig. 10) are sym¬ 
metric for both compact parallelogram shaped A- and 
B-type fee- and licp-pentamers as shown in Table III. 
Note that all concerted diffusion processes converts an 
A-type fcc-pentamer to a B-type hcp-pentamer (Fig. 10 
(c & d)) and vice versa. Similarly, all concerted diffusion 
processes converts a B-type fcc-pentamer to an A-type 
hcp-pentamer (Fig. 10 (a & b)) and vice versa. Hence if 
we start with an A-type fcc-pentamer, then only B-type 
hcp-pentamer is accessible through concerted processes 
whereas B-type fee- and A-type hep-pentamers are not 
accessible. Same is the case when we start either with 
B-type fcc-pentamer or A-type hcp-pentamer. This re¬ 
sults in anamolous diffusion for compact pentamer. In 
order to get random diffusion for compact pentamer, 
shape changing single-atom corner rounding processes 
are cruical. One such example of shape changing single 
atom corner rounding process is shown in Fig. 11 for an 
fee pentamer. To get an A-type fcc-pentamer (Fig. ll(b 
& c)), energy barrier for this single atom corner round¬ 
ing process is 0.136 eV with a reverse barrier of 0.396 
eV while to get a B-type fcc-pentamer (Fig. ll(b & a)), 
it is 0.063 eV with a reverse barrier of 0.315 eV. Simi¬ 
lar single atom corner rounding processes to get A- and 
B-type hep-pentamers from an hcp-pentamer were also 
found automatically during simulations. Our database 
for pentamer also contains various types of multi-atom 


processes found automatically during simulations. 




FIG. 10. Shows concerted diffusion processes for 5-atom com¬ 
pact clusters and their activation barriers 


TABLE III. Activation barriers (in eV) of concerted transla¬ 
tion processes for compact A- (and B-type) FCC pentamers 
and those for compact A- (and B-type) HCP pentamers. 

Directions fee hep 

f 0.288 (0.281) 0.286(0.279) 

2 0.290 (0.282) 0.287(0.281) 

3 0.290 (0.282) 0.287(0.281) 



FIG. 11. Interconvert process for 5-atom compact clusters 
and their activation barriers, change actual barriers 


5. Hexamer 

Compact B-type hcp-pentamer of Fig. 9(b) results in 
three compact liexamers, namely: A triangular B-type 
fee-hexamer (with all B-steps each of length 2 as shown in 
Fig. 12(b)) and two parallelogram shaped fcc-hexamers 
(symmetric) with same lengths of A and B-type steps 
(we will call them parallelogram fcc-hexamer) as shown 
in Fig. 12(c & d). From a compact A-type fcc-pentamer 
of Fig. 9(c), three compact liexamers, namely: A trian¬ 
gular A-type fcc-hexamer (with all A-steps each of length 



















FIG. 12. Compact hexamer shapes obtained when an atom is 
added to (a) a B-type fcc-pentamer at three different positions 
marked by b, c and d. (b) Compact triangular B-type fcc- 
hexamer with all B-steps. (c) Compact parallelogram shaped 
fcc-pentamer with same lengths of A- and B-steps, (d) Com¬ 
pact parallelogram shaped fcc-pentamer with same lengths 
of A- and B-steps. Note that compact parallelogram shaped 
hexamers shown in (c) and (d) are symmetric. 



FIG. 13. Activation barriers along direction 1 for the con- 
ceretd diffusion process of a compact parallelogram (a) fcc- 
hexamer (b) hcp-hexamer. 


TABLE IV. Activation barriers (eV) of concerted hexamer 
translations processes shown in Fig. 13 


Directions 

fee 

hep 

1 

0.306 

0.304 

2 

0.353 

0.351 

3 

0.244 

0.242 


2 as shown in Fig. 14(a)) and two parallelogram shaped 
fcc-hexamers (symmetric) with same lengths of A- and 13- 
type steps which are similar to those obtained from com¬ 
pact B-type fcc-pentamer. Similarly, each of the compact 
A- and B-type hcp-pentamers of Fig. 10(b & d) results in 
two compact triangular hcp-hexamers of A- (with all A- 
steps as shown in Fig. 14(d)) and B-type (with all B-steps 
as shown in Fig. 14(b)) respectively and two parallel¬ 
ogram shaped compact hcp-hexamers (symmetric) with 
equal lengths of A- and B-steps (as shown in Fig. 13(b)). 

We found that compact parallelogram fcc-hexamer 
with similar A- and B-steps (Fig. 13 (a)) is most stable 
whereas compact parallelogram hcp-hexamer with simi¬ 
lar A- and B-steps (Fig. 13 (b)) is 0.002 eV less stable. 
Also compact triangular A-type fcc-hexamer (Fig. 14 (a)) 
is 0.015 eV more stable than compact triangular B-type 
hcp-trangular hexamer (Fig. 14 (b)). whereas compact 
triangular A-type hcp-hexamer (Fig. 14 (d)) is 0.010 eV 
more stable than compact triangular B-type fcc-hexamer 
((Fig. 14 (c))).Fig. 13 shows concerted diffusion process 
for compact parallelogram fee- and hcp-hexamer along 
direction 1, activation barriers in all three directions are 
given in Tables. IV. As can be seen from Table. IV, the 
concerted diffusion barrier for compact parallelogram fcc- 
and hcp-hexamers of Fig. 13 are different in all three di¬ 
rections. Reported in Fig. 14 are the concerted diffu¬ 
sion processes for compact triangular A- and B-type fcc- 
and hcp-hexamer along direction 1. Since all three direc¬ 
tions for the compact triangular A- and B-type hexamers 
shown in Fig. 14 are symmetric, so their activation barri¬ 
ers for concerted diffusion are same in all three directions. 


6. Hept.amer 

On fcc(lll) surface heptamer has a compact closed- 
shell structure (having same A- and B-type steps of 
equal lengths) with each edge atom having at least three 
nearest neighbor bonds as shown in Fig. 15. There are 
two compact hexagonal shaped heptamers, one occupy¬ 
ing fee (called fcc-heptamer) and other occupying hep 
(called hep-heptamer) sites. Our SLKMC simulations 
found that both fee- and hep-heptamer diffuses primar¬ 
ily via concerted diffusion processes, which displace the 
cluster from fcc-to-hcp and vice versa, and their barriers 
are shown in Fig. 15. This can also be concluded from 
the effective activation barrier shown in table. 27, which 
is closer to the average of activation barriers shown in 
Fig. 15. Since the compact heptamer has a symmetric 
shape accordingly activation barriers in all the other di¬ 
rections are also the same as shown in Fig. 15. We note 
that fee heptamer is 0.002 eV more favorable than hep 
heptamer. Our database for hepatmer also contains var¬ 
ious single- and multi-atom processes as well. 

7. Octamer 

Adding another atom to the compact fcc-heptamer, 
we get two compact octamers, one with overall long B- 
type step egdes than A-type step edges (called B-type 
fcc-octamer) obtained by attaching an atom to one of 
the A-steps of fcc-heptamer (Fig. 16(a & b)) and an¬ 
other one with overall long A-type step egdes than B- 
type step edges (called A-type fcc-octamer) obtained by 
attaching an atom to one of the B-type step edges of 













9 


0.0 eV 






FIG. 14. Triangular shaped 6-atom island obtained by adding 
an atom to shorter edge of a 5-atom cluster, (a) all-fcc with 
B-type step edges (b) all-hep with A-type step edges (c) all- 
fcc with A-type step edges (d) all-hep with B-type step edges. 
Also shows activation barrier for concerted diffusion in direc¬ 
tion 1. 



FIG. 16. Compact geometries for fcc-octamers obtained by 
attaching an atom to (a) a compact fcc-heptamer at posi¬ 
tions marked as b and c with dotted circles, (b) fcc-octamer 
with overall long B-type step edges than A-type step edges 
(called B-type fcc-octamer). (c)fcc-octamer with overall long 
A-type step edges than B-type step edges (called A-type fcc- 
octamer). 



FIG. 15. Concerted diffusion processes and their activa¬ 
tion barriers in direction 1 for a (a) fcc-heptamer (b) hcp- 
heptamer. Note that due to symmetry, activation barriers 
along direction 2 and 3 are same as those along direction 1 
for the two types of heptamer. 


fcc-heptamer (Fig. 16(a & c)). Similarly two compact B- 
and A-type hcp-octamer can be obtained by attaching 
another atom to any of the A- and B-type step edges of 
the compact hcp-heptamer of Fig. 15 are shown in Fig. 17 
(b & d). Compact A-type fcc-octamer is 0.008 eV more 
stable than compact B-type hcp-octamer whereas com¬ 
pact A-type hcp-octamer is 0.002 eV more stable than 
compact B-type hcp-octamer. 

Compact octamers diffuse via concerted diffusion pro¬ 
cesses which displaces the island from all-fcc sites to all- 
hcp sites and vice versa as shown in Fig 17 (a - d). In 
general, activation barriers for compact A- and B-type 


fcc-octamers are different. Same is the case for A- and 
B-type hcp-octamers. It can be seen in the Fig. 17 (a - d) 
that concerted diffusion process convert an A-type fcc- 
octamer into a B-type hcp-octamer and vice versa. Ta¬ 
ble. V gives activation barriers for concerted diffusion 
processes in all 3 directions for both compact fee- and 
hcp-octamers (see Fig. 17). We note that concerted dif¬ 
fusion processes in the directions 2 and 3 are symmetric 
and have lowest energy barrier whereas for direction 1 it 
is highest. 

Although compact octamers diffuse primarily via con¬ 
certed diffusion processes, we found in our simulations 
that both multi-atom (including shearing, reptation, 
rotation etc) and single-atom (including kink attach¬ 
ment/detachment, corner rounding etc.) processes are 
also relatively common. Accordingly we will discuss 
multi-atom processes particular to this island size here, 
while single atom processes, which are common to all 
island sizes, will be discussed in detail later. 

Interesting multi-atom processes includes shearing and 
reptation as shown in Figs. 18 & 19 are also found in our 
simulations. In the case of shearing processes, part of 
the island (more that one atom) moves from fee (hep) 
to the nearest fee (hep) sites if the island is initially on 
fee (hep) sites. An example of trimer shearing process 
along with their activation barriers is shown in Fig. 18 
(a) & (b) while Fig. 18 (c) & (d) shows a dimer shearing 
process and its activation barrier. 
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TABLE V. Activation barriers (eV) for concerted diffusion 
processes for compact A-type (B-type) fcc-octamers and com¬ 
pact B-type (A-type) hcp-octamers as shown in Fig. 17 (a - 

d). 


Directions 

fee 

hep 

1 

0.472 (0.467) 

0.463 (0.469) 

2 

0.418 (0.411) 

0.409 (0.413) 

3 

0.418 (0.411) 

0.409 (0.413) 

0.0 eV 

(b) 

0.009 eV 




FIG. 17. Concerted diffusion processes along direction 2 for a 
compact (a) A-type FCC octamer (b) B-type HCP octamer 
(d) B-type FCC octamer and (d) A-type HCP octamer. 


Another important multi-atom process is reptation 
through which an entire island diffuses from all-fcc (all- 
hcp) to all-hep (all-fcc) sites in a two step shearing pro¬ 
cess. In the first step part of the island moves from fee 
(hep) to nearest hep (fee) sites (see Fig.l9(a & b)) and in 
this intermediate state, part of the island is on fee sites 
and part of it is on hep sites. In the next step, part of 
the island which is on fee (hep) sites moves to hep (fee) 
sites (see Fig.l9(a & b)) thus completing the reptation 
mechanism. 


8. Nonamer 

By addition of another atom to either of the com¬ 
pact A- or B-type fcc-octanrer results in a single dia¬ 
mond shaped compact nonamer as shown in Fig. 20(a 
- b) for the case of B-type fcc-octamer to a compact 
diamond shaped fcc-nonamer (having equal A- and B- 
step lengths). Similarly, compact A- or B-type hep- 
octanrer results in a single diamond shaped compact hep- 
nonamer. Two other almost compact diamond shaped 
configurations of a nonamer for each fee- and hep- non- 



FIG. 18. Examples of shearing diffusion processes along with 
their activation barriers for an 8-atom islands, (a & b) Trimer 
shearing (c & d) Dimer shearing. 



FIG. 19. Shows different processes (or step) involved in 8- 
atom reptation diffusion mechanism 




FIG. 20. Compact diamond shaped nonamer obtained by 
(a) Adding an atom to the compact B-type fcc-octamer (b) 
Resultant compact diamond shaped fcc-nonamer. 





















11 


0.0 eV -0.002 eV 



FIG. 21. Concerted diffusion processes and their activation 
barriers for 9-atoms compact islands. 

TABLE VI. Activation barriers (eV) of concerted translations 
processes in all 3 directions for the 9-atom cluster shown in 
figs. 21 (a & b, c & d) 


Directions 

fee 

licp 

1 

(0.518, 0.414) 

(0.516, 0.411) 

2 

(0.474, 0.526) 

(0.472, 0.523) 

3 

(0.509, 0.414) 

(0.507, 0.411) 


amer with a kink in A- or B-step (See Fig. 21 (a - b)) 
are also oberved in out simulations. We find that fcc- 
nonamer with a kink in A-step is most stable configura¬ 
tion whereas fcc-nonmaer with a kink in B-step (Fig. 21 
(d)) is 0.002 eV higher in energy. Similarly hcp-nonamers 
with a kink in A-step and B-step (Fig. 21 (a)) are less 
stable than most stable nonamer by 0.004 eV and 0.005 
eV respectively. Least stable configuration are the sym¬ 
metric fee- and hcp-nonamers (Fig. 21 (e & f)) by 0.007 
eV and 0.010 eV respectively. 

Most frequently observed processes for a nonamer are 
single-atom processes like edge diffusion processes either 
along A- or B-type step edges or corner rounding and 
kink detachment processes (K3 —» A2(B2), see Fig. 25). 
In contrast, concerted diffusion processes although picked 
less often but does contribute most to the nonamer diffu¬ 
sion, since they produce large displacement in the center 
of mass, and as usual with these types of processes a 
compact island converts from an all-fcc to an all-hep is¬ 
land. Although single-atom process individually do not 
produce big displacement to the center of mass of an is¬ 
land, in case of a nonamer they are the most frequently 
observed processes and do contribute to smaller extent to 
the island diffusion. This can be concluded by comparing 
effective activation barrier for this island size in table. IX, 
and activation barriers for single atom processes in ta¬ 
ble. VIII and concerted processes in table. VI. Effective 
activation barrier is slightly larger than the maximum 
activation barrier for most frequently observed concerted 
processes. Some of the most frequently observed compact 


0.0 eV 0.0 eV 




FIG. 22. Dimer shearing process and their activation barriers 
for compact 9-atom island. 


nonamer shapes, corresponding concerted diffusion pro¬ 
cesses and their activation barriers are shown in Fig. 21. 
Note that energy barriers for fee (hep) clusters shown in 
figs. 21(a - d) in directions 1 and 3 have same energy 
barrier. Table. VI shows the activation barriers for con¬ 
certed processes in all 3 directions for nonamers shown in 
Fig. 21. Furthermore for these 9-atom islands, concerted 
diffusion in the direction 1 & 2 are the most frequently 
observed concerted processes during the island diffusion. 

Fig. 22 shows frequently observed multi-atom process 
and its activation barrier during nonamer diffusion. This 
is again is a shearing of dimer along A-type step-edge 
similar to what was observed previously with other is¬ 
land sizes. With this shearing process atoms diffuse from 
one fee (hep) to another fee (hep) sites and appears as 
though a dimer is diffusion along an edge of compact 
7-atom cluster. Activation barrier for this shearing pro¬ 
cess is lower than a single-atom diffusion along an edge 
as well as a corner rounding process. We note that this 
dimer forms an A-type step edge to the 9-atom island. 
Although not shown here, we have also observed repta- 
tion processes 3 ' as well and activation barriers for the 
intermediate processes were ever lower than single-atom 
diffusion processes, but former happens only when the 
shape of the island is non-compact. 


9. Decamer 

There are two possible compact geometries, one for 
each fee- and hcp-decamer as shown in Figs. 23(a & b) 
with fcc-decamer being 0.003 eV more stable than com¬ 
pact hcp-decamer. In the case of a decamer we have 
observed a variety of single-atom, multi-atom and con¬ 
certed diffusion processes as well. Most frequent shaped 
observed in our simulations are the compact shapes as 
shown in Fig. 23 which has same number of A- and B- 
type step edges of same length. For this shape most 
frequently observed process was concerted diffusion pro¬ 
cesses as shown in Fig. 23 along with their activation 
barriers and as usually changes all-fcc cluster into a all- 
hep cluster and vice versa. Table. VII shows activation 
barriers for the concerted process shown in Fig. 23 in all 
3 directions both for fee and hep clusters. 

For a non-compact 10-atom cluster we have also ob- 
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FIG. 23. Concerted diffusion processes along direction 1 for 
(a) fee decamer (b) hep decamer. Values in brackets are the 
energy barriers obtained from DFT calculations. 


TABLE VII. Activation barriers (eV) of concerted cluster 
translations processes in all 3 directions for the compact fee 
and hep decamers of Figs. 23. 


Directions 

fee 

hep 

1 

0.580 

0.576 

2 

0.559 

0.555 

3 

0.559 

0.555 


served multi-atom processes like shearing and reptations. 
Dimer shearing process along A-type step-edge similar to 
what has been observed for small clusters is the most fre¬ 
quently observed shearing process. Fig. 24 shows diffu¬ 
sion processes involved in a reptation process along with 
their activation barriers 


10. Single-atom processes 

Various types of single-atom processes like edge dif¬ 
fusion, corner rounding, kink attachment, kink detach¬ 
ment etc were automatically identified and their activa¬ 
tion barriers were found during our simulations. Some 
of these single-atom processes are shown in Fig. 24(a & 
b) for an hep island with their activation barriers and 
those for their fee analogs are given in Table VIII. In 



FIG. 24. A step in reptation process for a decamer along with 
its activation barrier, (a) five atom island moving from fee to 
hep sites, (b) its reverse. 


each single-atom process, an atom on an fee site moves 
to a nearest-neighbor vacant fee site, while an atom on 
hep site moves to a nearest-neighbor vacant hep site. 
These single-atom processes are also referred as long- 
jumps. We also found, specially for small islands of size 
1-6, single-atom processes, called short-jumps, in which 
an atom on an fcc(hcp) site moves to nearest-neighbor 
hcp(fcc) site. These short-jump single-atom processes are 
dominant diffusion mechanism for hetero-epitaxial sys¬ 
tem with large lattice mismatch. Note that the activa¬ 
tion barriers for single-atom processes depend not only 
on whether the atom is part of an fee island or an hep 
island but also on the neighborhood of the diffusing atom 
(like atom diffusing on an A-type or a B-type step edge, 
corner, kink etc). 

To classify single-atom processes in Table VIII we have 
used the notation X ni U -A Y n/ V. where X or Y = A 
(for an A-type step-edge) or B (for a B-type step-edge) 
or K (for kink) or C (for corner) or M (for monomer); 
n i = the number of nearest-neighbors of the diffusing 
atom before the process; n f = the number of that atom’s 
nearest neighbors after the process. U or V = A or B 
(for corner or kink processes) or null (for all other other 
process types). 

For example, process 1, B 2 —> B 2 , is a single-atom 
B-step edge process in which the diffusing atom has 2 
nearest-neighbors before and after the process. Process 
3, CiB -A B 2 , is a corner rounding process towards a 
B-step, the diffusing atom starting on the corner of a 
B-step with one nearest-neighbor and ending up on the 
B-step with two nearest-neighbors. In process 10, C 2 A 
—> CiB, the diffusing atom begins on the corner of an 
A-step having two nearest-neighbors and ends up on the 
corner of a B-step with only one nearest-neighbor. 


IV. DIFFUSION COEFFICIENTS AND 
EFFECTIVE ENERGY BARRIERS 


We start our simulations with an empty database. Ev¬ 
ery time a new configuration (or neighborhood) is found, 
SLKMC-II finds all possible processes and their acti¬ 
vation barriers for it on-the-fly and stores them in a 
database. Since the processes and their activation bar¬ 
riers are dependent on island size, each one of them 
has separate unique database. We carry out 10 7 KMC 
steps for each island size at temperatures 300K, 400K, 
500K, 600K and 700K. The diffusion coefficient of an 
adatom island is calculated by D = linp^oo ([RcomW 
- Rcom( 0)] 2 )/2dt, where D is the diffusion coefficient, 
RcOM(t) is the position of the center of mass of the is¬ 
land at time t, and d is the dimensionality of the system 
which in our case is 2. Diffusion coefficients thus obtained 
for island sizes 1 — 10 are summarized in Table IX. Ef¬ 
fective energy barriers for each island size are extracted 
from their respective Arrhenius plots as shown in Fig. 26 
while Fig. 27 shows plot of effective energy barriers as a 
function of island size. 
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TABLE VIII. Activation barriers (eV) of single-atom pro¬ 
cesses for both fee and hep islands. The index numbers refer 
to the types of processes illustrated in Fig. 25. See text for 
explanation of the notation used to classify the process types. 


Index no. 

Process type 

fee 

hep 

1 

B 2 —>• B 2 

0.349 

0.347 

2 

B 2 —y M 

0.555 

0.555 

3 

CiB -> B 2 

0.138 

0.071 

4 

CiB -> M 

0.319 

0.319 

5 

CiB -» A 2 

0.72 

0.138 

6 

CiB -> CiB 

- 

- 

7 

CiB -> M 

- 

- 

8 

A 2 —y M 

0.537 

0.537 

9 

A 2 —y A 2 

0.266 

0.265 

10 

C 2 A ->• CiB 

0.385 

0.308 

11 

C 2 A ->• M 

0.564 

0.545 

12 

B 2 —y K 3 B 

0.322 

0.322 

13 

K 3 A —> a 2 

0.471 

0.470 

14 

a 2 k 3 a 

0.252 

0.252 

15 

K 3 B — y B 2 

0.535 

0.535 

16 

C 3 B -> CiA 

0.579 

- 

17 

C 3 B -»■ M 

0.753 

- 

18 

C 3 B -> CiB 

0.516 

- 



(b) 



FIG. 25. (Color online) Single-atom processes for an hep is¬ 
land. The index numbers designate the processes described 
in Table. VIII, which gives the activation barriers for each. 



FIG. 26. Arrhenious plots for 1—10 atom islands. 



FIG. 27. Effective energy barriers of 1 — 10 atom islands as a 
function of island size. 

V. DISCUSSION 

To summarize, we have performed a systematic study 
of the diffusion of small Ag (1-10 atom) islands on 
Ag(lll) surface, using newly developed self learning 
KMC (SLKMC-II) simulations in which the system is al¬ 
lowed to evolve through mechanisms of its choice with the 
usage of a self-generated database of single-atom, multi¬ 
atom and concerted diffusion processes. Specifically, dif¬ 
fusion processes like concerted processes and multi-atom 
processes involving fcc-fcc, fcc-hcp and hep-hep jumps 
are automatically found during the simulations. We find 
that these small-sized islands diffuse primarily through 
concerted motion with a small contribution from single 
atom processes, even though for certain cases the fre¬ 
quency of single atom processes is large because of lower 
activation energies. Furthermore multi-atom processes 
including shearing and reptation, which are activated at 
higher temperature, contributes rarely to the diffusion of 
island. Interestingly enough, energy barriers for repta¬ 
tion processes are small but they are only occurred when 
island shape is non-compact which is rare for the tern- 
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TABLE IX. Diffusion coefficients ((A°) 2 /sec) at various temperatures and effective energy barriers for Ag islands. 


Island Size 

300K 

400K 

500K 

600K 

700K 

Eeff (eV) 

1 

2.11 

X 

10 11 

3.74 

X 

10 11 

6.62 

X 

10 11 

5.00 

X 

10 11 

7.79 

X 

10 11 

0.059 

2 

6.33 

X 

io 10 

1.68 

X 

io 11 

2.76 

X 

io 11 

4.13 

X 

io 11 

5.29 

X 

io 11 

0.096 

3 

2.94 

X 

10° 9 

1.62 

X 

io 10 

4.72 

X 

io 10 

9.51 

X 

io 10 

1.49 

X 

io 11 

0.179 

4 

2.20 

X 

1 0 °9 

1.57 

X 

io 10 

4.48 

X 

io 10 

1.14 

X 

io 11 

1.86 

X 

io 11 

0.202 

5 

1.08 

X 

10°8 

1.72 

X 

io 09 

8.73 

X 

io" 

2.83 

X 

io 10 

6.01 

X 

io 10 

0.287 

6 

4.50 

X 

10° 7 

7.87 

X 

10°8 

4.36 

X 

io" 

1.45 

X 

io 10 

3.29 

X 

io 10 

0.298 

7 

1.14 

X 

1 0 °6 

6.14 

X 

IO 07 

6.90 

X 

l 0 O8 

3.07 

X 

io 09 

1.01 

X 

io 10 

0.410 

8 

9.55 

X 

IQ 04 

4.99 

X 

10°6 

5.64 

X 

IO 07 

2.81 

X 

io 08 

1.00 

X 

i 0 °9 

0.413 

9 

3.48 

X 

IQ 04 

2.07 

X 

10°6 

2.51 

X 

IO 07 

1.90 

X 

io 08 

7.12 

X 

io°s 

0.425 

10 

2.15 

X 

io 02 

7.49 

X 

IO 04 

2.51 

X 

IO 06 

2.57 

X 

IO 07 

1.29 

X 

10°8 

0.501 



FIG. 28. Distribution of single-atom, multi-atom, concerted 
and total processes for 1-10 atom islands accumulated in the 
database during SLKMC simulations. 


perature range we study here. We present in Fig. 28 the 


total number of processes found by our simulations for 
different island sizes. As can be seen from it, number 
of processes increase exponentially, specially beyond 4- 
atom island, suggesting the need for an automatic way 
of finding all the possible processes during simulations. 
Also by the inclusion of fee and hep sites in the pattern 
recognition scheme allows us to get all possible concerted 
processes which involve motion from fcc-hcp sites. Fig. 28 
also shown single atom and multi atom processes. Allow¬ 
ing the system the possibility of evolving in time through 
all types of processes of its choice, we are able to estab¬ 
lish the relative significance of various types of atomistic 
processes through considerations of the kinetics and not 
just the energetics and/or the thermodynamics, as is of¬ 
ten done. For small Ag islands on Ag(lll), we find the 
effective barriers for diffusion to scale with island size. 
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